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Using first-principles calculations, the effect of magnetic point defects (vacancy and adatom) is in- 
vestigated in zigzag graphene nanoribbons. The structural, electronic, and spin-transport properties 
are studied. While pristine ribbons display anti-parallel spin states at their edges, the defects are 
found to perturb this coupling. The introduction of a vacancy drastically reduces the energy differ- 
ence between parallel and anti-parallel spin orientations, though the latter is still favored. Moreover, 
the local magnetic moment of the defect is screened by the edges so that the total magnetic moment 
is quite small. In contrast, when an adatom is introduced, the parallel spin orientation is preferred 
and the local magnetic moment of the defect adds up to the contributions of the edges. Furthermore, 
a spin-polarized transmission is observed at the Fermi energy, suggesting the use of such a defective 
graphene nanoribbon as spin- valve device. 

PACS numbers: 61.46.-w, 71.15.Mb, 73.63.-b, 75.75.+a 



Recent progresses in the preparation of contacted 
graphene mono-layer have enabled the synthesis of well- 
controlled nanometer-sized systems with open bound- 
aries [H 12] . Both one-dimensional graphene ribbons and 
zero-dimensional graphene dots have been achieved ei- 
ther by cutting exfoliated graphene layers [3l H] or by 
patterning epitaxially grown graphenes [5]. Chemical 
methods have also been designed to achieve a solution- 
phase derivation of graphene ribbons with ultrasmooth 
edges [6 . These advances are at the origin of the rais- 
ing interest in open boundary systems where the pres- 
ence and shape of the edges influence drastically the tt- 
electronic structure [7], opening the road towards new 
physical phenomena and technological applications. 

As far as ribbons are concerned, the geometry of sp'^- 
bonded networks implies two possible cutting directions 
called "zigzag" and "armchair" according to the shape 
of the created edge. Due to topological reasons, zigzag- 
shaped edges give rise to peculiar extended electronic 
states which decay exponentially inside the graphene 
sheet [H E]. These edge-states, which are not reported 
along the armchair-shaped edges, come with a twofold 
degenerate flat band at the Fermi energy over one third 
of the Brillouin zone and have been suggested to be 
ferromagnetically ordered along the edge [10]. Recent 
DFT calculations on GNRs with zigzag-shaped edges 
(zGNRs), have suggested that due to the finite width 
of the ribbon, the interaction between the edges favors 
an opposite spin orientation at these edges ^1 |Tl] [12] . 
These predictions are in agreement with the second Lieb's 
theorem [13] about the total spin in bipartite lattices. 

Besides the possible magnetization of the edges, defect- 
induced magnetism is also expected in GNRs. In- 
deed, while many common topological defects (5/7 pairs. 
Stone- Wales [14 ,.••) do preserve the saturation of the 
carbon atoms and are reported not to cause any lo- 



cal magnetization [15], processes where dangling bonds 
are created, may possibly induce a local magnetic mo- 
ment [ini [TTl [18]. For instance, carbon atoms can be 
removed from their lattice position, creating pairs of va- 
cancies and adatoms. Such Frenkel defects are frequently 
present at thermal equilibrium [19,20^ in sp^-bonded car- 
bon systems. Depending on their positions, such mag- 
netic point defects might affect differently the magnetic 
structure of GNRs. When located at the edge, these 
defects affect mainly the ferromagnetic order along the 
edge [21]. However, when located inside the ribbon, their 
effect on the magnetic structure is still an open issue. 

In the present letter, magnetic point defects located 
in the vicinity of the ribbon axis have been used to 
tailor both the electronic and spin-transport properties 
of hydrogen passivated zGNRs. More specifically, the 
introduction of vacancy /adatom in these systems de- 
creases the energy difference between the parallel and 
anti-parallel spin orientations of the edges. The mag- 
netic moment of the vacancy defect is found to be dras- 
tically reduced in the defective ribbon compared to the 
ideal one. In the presence of carbon adatoms, a paral- 
lel spin orientation at the edges is favored, giving rise to 
a spin-polarized electronic transmission function around 
the Fermi energy. 

Electronic irradiations can induce such magnetic point 
defects in graphene samples and allow to tune their num- 
ber and their concentration [22]. Therefore, based on 
our theoretical predictions, one can hope to use defective 
zGNRs as spin-based nano-electronic devices. 

Spin Polarized Density Functional Theory [23l [24] 
within the generalized-gradient approximation is used 
to investigate a hydrogen passivated 8-zGNR (width 
~ ISA). Periodic boundary conditions are used with fixed 
lateral dimensions which insure 25 A of vacuum between 
the GNRs in neighboring cells. When vacancies or car- 
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bon adatoms are introduced, a 6x1x1 supercell is con- 
sidered, leading to an nearest-neighboring defect distance 
of ~ 15 A and a defect concentration of ~ 1%. Note 
that, for the computation of the transmission functions, 
a 22x1x1 superceh is considered in order to obtain a 
good screening of the perturbed Hartree potential due 
to the defect. In order to deal with the large number of 
atoms in the supercells, numerical atomic orbital basis 
set [25] are used to expand the wave- functions in con- 
junction with norm-conserving pseudo-potentials. The 
energy levels are populated using a Fermi-Dirac distribu- 
tion with an electronic temperature of 250K. The inte- 
gration over the ID Brillouin zone is replaced by a sum- 
mation over a regular grid of 40 A:-points [26 . The geom- 
etry is fully relaxed until the forces on each atom, and 
on the unit cell [27] are less than 0.01 eV/A and 0.05 
eV/A respectively. The small spin-orbit coupling of car- 
bon atoms is neglected and independent collinear spin 
orientations, labeled a and are considered. Within 
such framework, the electronic ground state of the ideal 
8-zGNR reveals anti-parallel (||) spin orientations be- 
tween the edges leading to a semiconducting band struc- 
ture (0.5 eV band gap) with full spin degeneracy. On 
the other hand, the magnetic configuration with parallel 
(It) spin orientations between the edges is found to be 
metastable (11 meV/edge-atom higher in energy). This 
configuration displays a metallic behavior, as the tt* and 
7r/3 bands are crossing at the Fermi energy, and presents a 
total magnetic moment of 0.51 fiB per edge atom. These 
results are in good agreement with the previous calcula- 
tions [101 ng. 

By analogy with graphene, the removal of a carbon 
atom (vacancy creation) induces a Jahn- Teller distortion 
of the honeycomb structure. Two of the unsaturated 
carbon atoms come closer one to the other and form a 
weak covalent bond, inducing a pentagon like rearrange- 
ment. The third unsaturated carbon atom moves out of 
the plane [161 EH]- Consequently, the initial D^^ sym- 
metry of the hexagonal network is lost in favor to the 
Cs symmetry [29]. In zGNRs, such atomic rearrange- 
ment leads to two possible orientations of the vacancy, 
called "perpendicular" and "tilted" according to the rel- 
ative orientation of the mirror plane of the vacancy with 
respect to the ribbon axis. 

As expected due to the presence of dangling bonds, the 
electronic density appears to be locally spin-polarized on 
the vacancy. The interaction between the magnetic mo- 
ment of the defect and the spin ordered edge states results 
in four locally stable configurations for both vacancy ori- 
entations (see Fig.jl]). In the following, the various mag- 
netic configurations will be labeled Xy, where X refers 
to the parallel (||) or anti-parallel (ti) spin orientations 
between the ribbon edges, while Y refers to the parallel 
(P) or anti-parallel {AP) spin orientations between the 
defect and the nearest edge. 

The energies related to the eight possible magnetic con- 




FIG. 1: [Color online] Spin polarized electronic densities 
{pa — P(3) isosurfaces for the four magnetic configurations as- 
sociated to the tilted vacancy in zGNR. Green and yellow 
surfaces correspond to an excess of spin a and spin [3 elec- 
trons, respectively. 



Tilted vacancy Perpendicular vacancy 
(eV) Mtot (pb) AE (eV) Mtot (pb) 



Hap 


0.000 


0.19 


0.110 


1.98 


Up 


0.004 


1.68 


0.142 


0.06 


TTap 


0.021 


2.03 


0.154 


2.10 


TTp 


0.023 


4.06 


0.170 


4.52 



TABLE I: Energies related to the ground state (TIap) and 
total magnetic moments associated to the eight possible mag- 
netic configurations of the vacancy. 



figurations of the vacancy are detailed in Table |l| The 
tilted vacancy is found to be the most stable orientation, 
probably because of the formation of a stronger C — C 
bond during the atomic reconstruction process. Con- 
sequently, the ribbon width displays a small shrinking 
(1.1%) around the vacancy, while a stretching is observed 
for the perpendicular vacancy orientation. Both vacancy 
orientations have the same ][ ground state as the pristine 
GNR. However, the energetics of the inter-edge coupling 
is modified, as shown in Table [Tj Indeed, while the energy 
difference between the parallel (||) and anti-parallel (ti) 
spin orientations at the edges is 66 meV in the absence 
of defect, this value is reduced to 21 meV (44 meV) in 
the presence of a tilted (perpendicular) vacancy. At last, 
an accurate convergence study with respect to the size of 
the supercell allows us to predict a formation energy of 
15.6 eV (which is approximately twice the one reported 
for the single vacancy in the graphene sheet [iF, T8^) and 
a weak magnetic moment of 0.15 fiB for the isolated va- 
cancy in the tilted Hap configuration. 

At equilibrium, the carbon adatom is in a bridge-like 
position between two in-plane carbon atoms. This geom- 
etry is similar to the one reported for carbon atoms on 
graphene [H HH IH [30l [31] . Again, in zGNRs, two ad- 
sorption sites are possible: labeled perpendicular or tilted 
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FIG. 2: [Color online] Spin polarized electronic densities {pa — 
pjs) isosurfaces for the four magnetic configurations associated 
to the carbon adatom in the tilted orientation. Green and 
yellow surfaces correspond to an excess of spin a and spin [3 
electrons, respectively. 





Tilted adatom 


Perpendicular adatom 




(eV) Mtot (pb) 


(eV) 


Mtot (pb) 


TTp 


0.000 3.73 


0.060 


4.00 


TTap 


0.005 2.58 


0.089 


2.81 


Tip 


0.008 0.11 


0.094 


0.72 


TIap 


0.032 0.67 


0.094 


0.72 



TABLE II: Energies related to the ground state (TTp) and 
total magnetic moments associated to the eight possible mag- 
netic configurations of the carbon adatom. 



depending on the orientation of the underlying carbon- 
carbon bond with respect to the ribbon axis. In both 
cases, the carbon adatom displays a local magnetic mo- 
ment. In analogy with the vacancy case, the various mag- 
netic configurations are labeled using the aforementioned 
Xy syntax (see Fig. [2|, and their corresponding energies 
are detailed in Tablellll 

The tilted C — C bond is found to be the preferred host 
for carbon adsorption as it is associated with a smaller 
in-plane C — C bond and a larger distance between the 
adatom and the ribbon. No significant shrinking is re- 
ported here. Similar to vacancies, carbon adatoms mod- 
ify the energetics of the inter-edge coupling. However, 
this effect is enhanced in the presence of carbon adatoms 
and the TT ground state is different from the one the 
pristine GNRs. Finally, an adsorption energy of 1.2 eV 
and a magnetic moment of 0.48 fiB are computed for the 
isolated carbon adatom in the TTp configuration. These 
values are quite similar to the 1.40 eV and 0.45//^ re- 
ported for the adsorption of carbon on graphene [16]. 

While the resulting ground states are different in the 
presence of a vacancy (Ti) or a carbon adatom (TT), both 
defects induce the same trend on the inter-edge coupling. 
Indeed, the energy difference between the Ti and the TT 
configurations in the presence of defects is found to be 



reduced well below the 66 meV predicted for the pristine 
ribbon. This result is a direct consequence of the defect- 
induced perturbation of the spin-polarized density. On 
one hand, by breaking the bipartite character of the car- 
bon network, defects partially destroy the enhanced ex- 
change splitting responsible of the stabilization of the Ti 
configuration [10]. On the other hand, as can be seen 
from Figs. [T] and [2] the spin polarization is larger in the 
center of the ribbon for the Ti configuration than for 
the TT orie. As a result, the density rearrangement is 
more limited around the defects in the former, which is 
confirmed by the projected density of states (not shown 
here). In summary, the introduction of defects induces a 
stabilization of the TT configuration in contrast with the 
pristine case. 

Since the introduction of magnetic point defects in 
zGNRs favors a specific spin configuration of the edges, 
a major impact is also expected on the transport proper- 
ties of these ID systems [32]. This will be investigated in 
the remainder of this letter. Our ab initio calculations of 
the electronic transmission functions, reported in Fig. [3j 
are performed within the Landauer approach. In order 
to simulate open boundary conditions, self-energies asso- 
ciated with the leads are included in the self-consistent 
calculation of the potential [33 . The supercell contain- 
ing the point defect is connected to two leads consisting 
of a few unit cells of ideal zGNR [Fig. [sj^a)]. Note that, 
as indicated above, a larger supercell has been used in 
order to ensure an accurate alignment of the electronic 
levels of the lead with those of the supercell. 

As mentioned before, the Ti spin configuration of the 
pristine zGNR is preserved when an isolated vacancy is 
introduced. The defective zGNR remains semiconduct- 
ing and its electronic transmission function displays a gap 
of 0.5 eV around the Fermi energy. On the contrary, in 
the presence of a carbon adatom, the TT spin configura- 
tion is favored and the zGNR becomes metallic, inducing 
a non-zero electronic transmission function at the Fermi 
energy. 

The main impact of the magnetic point defects on the 
transport properties is a global reduction of the trans- 
mission associated with the tt and tt* electrons. This can 
be related to a decrease of the transmission probability 
of some TT-TT* conduction eigenchannels compared to the 
pristine zGNR. For the vacancy [Fig.[3];b)], this effect ap- 
pears essentially for energies ranging from -0.80 to -0.15 
eV (tt channels) and from 0.4 to 0.5 eV (tt* channels), 
inducing a slight breaking of the spin degeneracy. For 
the carbon adatom, a similar reduction of the transmis- 
sion is found. But, defect-induced drops also appear at 
-0.3 eV and 0.02 eV. At these energies, states localized 
on the defect are present, which can only mix with the 
spin a conduction channel. Consequently, its transmis- 
sion probability goes to zero due to the alteration of its 
TT* character and the spin-degeneracy of the electronic 
transmission function is lifted up just above the Fermi 
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(E-Ep) [eV] 



FIG. 3: [Color online] (a) Schematic representation of the 
model structure used to predict the electronic transmission 
functions of the zGNR in the presence of (b) an isolated va- 
cancy in the TIap configuration, and (c) an isolated carbon 
adatom in the Hp configuration. The spin a and spin [3 
components are appear in blue and red, respectively, while 
the dotted line indicates the number of tt channels. 



energy. 

In conclusion, the || spin configuration of zGNRs 
tends to be favored by the introduction of magnetic 
point-defects in the ID system, as suggested by the low- 
energy spin-polarized density. The major impact of these 
point defects on the transport properties of zGNRs has 
been investigated. The predicted spin-polarized trans- 
mission function underlines the possibility to use adatom- 
doped graphene nanoribbons as valuable nano-devices in 
future spin-based electronics. 
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